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ABSTRACT 

Cosmological fireball models of 7-ray bursts (GRBs) predict delayed emission, 
"after-glow," at longer wavelengths. We present several new results regarding the 
model predictions, and show that X-ray to optical observations of GRB970228 and 
GRB970402 are naturally explained by the model: The scaling of flux with time 
and frequency agrees with model predictions and requires a power law distribution 
of shock accelerated electrons dlog A''e/(ilog7e = 2.3 it 0.1 (implying, and consistent 
with the observed, decline of flux observed at a given frequency); The absolute 
flux value agrees with that inferred through the model from observed 7-ray fluence. 
The future after-glow emission of these bursts is predicted. The observations indicate 
that the ratio of magnetic field to equipartition value and the fraction of dissipated 
kinetic energy carried by electrons are not much smaller than 1. More frequent 
observations at a single wavelength, or a wide spectrum at a single time, would put 
strong constraints on these parameters. We show that inverse-Compton scattering 
suppresses X-ray/optical emission at delays t < tic = 10(^e/0.3)^hr. Observations 
therefore imply < 0.3. The strong dependence of tjc on implies that tic ™ay vary 
widely from burst to burst, and that frequent X-ray /optical observations at t ~ tic 
would determine (,e- For ~ 0.2, inverse-Compton emission dominates during the first 
2hr, producing > IGeV photons and providing a natural explanation to the delayed 
GeV emission observed in several strong bursts. 



Subject headings: gamma rays: bursts 



Introduction 



Recent observations of 7-ray bursts (GRBs) suggest that they originate from cosmological 
sources ( Meegan et al. 1992 , paczyhski 1992 ). General phenomenological considerations indicate 
that the bursts are produced by the dissipation of the kinetic energy of a relativistic expanding 



fireball (see Piran 1994, Meszaros 1995 for reviews). After producing the main GRB, the expanding 



cooling fireball is expected to produce delayed emission at longer wavelengths ( Paczyhski fc 



Rhoads 199S, Meszaros & Rees 1997). Such optical and X-ray after-glow has possibly been 



detected for GRB970228 and GRB970402: Variable X-ray sources have been detected within the 
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3' GRB error boxes, fading on a time scale of days ( posta et al. 1997a ,b, Feroci et al. 1997 , Pirc| 



et al. 1997); A variable optical source, fading on a similar time scale, was detected within the 50" 



error box of the X-ray source associated with GRB970228 ( Proof et al. 1997] , |Sahu et al. T997| ). 



The GRB fireball model depends on several parameters: E, the total fireball energy; 7]~^, 
the fraction of energy initially carried by (proton) rest mass; T, the energy emission duration; 
the fraction of dissipated kinetic energy carried by electrons; (,b, the ratio of magnetic field to 
equipartition field; p, the spectral index of shock accelerated electrons {dNe/d'je oc 7~^). In sec. 2 
below we derive the after-glow predictions of the fireball model. Our results follow closely those of 
Meszaros & Rees (1997), the main new results being: (i) After-glow due to interaction with the 
inter-stellar medium (ISM) is determined by E, ^e, (,B, and p only (in particular, it is independent 
of 7] and of whether the energy is emitted impulsively or as a wind): The flux peaks at a frequency 
Um = Ci£;i/2^2^Bi-3/2, and for u > Um = C2^CBn^'^^(i^/i^m)"", where a = (p - l)/2 and n 
is the ISM density; (ii) Although p is expected to be similar for the GRB and the after-glow, 
a during the after-glow is expected to be smaller by 1/2 due to increase in cooling time. For 
the GRB a ~ 1, implying a ~ 0.5 for the after-glow. This also implies a decline in flux at a 
fixed wavelength much slower than the decline obtained assuming similar a for the GRB 

and after- glow; (iii) Although the delayed X-ray/ optical emission must result from synchrotron 
emission of accelerated ISM electrons, inverse-Compton emission may dominate during the first 
hours-days, suppressing the X-ray/optical emission. 

Observations are compared with model predictions in sec. 3. In sec. 4 we summarize the 
comparison results, discuss the observational constraints on the model and point out future 
observations that may further constrain model parameters, give the model predictions for future 
GRB970228 and GRB970402 after-glow, and discuss the sensitivity of the results to underlying 
assumptions. 



2. Model and Parameters 

Whatever the sources of GRBs are, observations suggest the following scenario for the 
emission of the observed 7-rays. A compact, rg ~ lO^cm, source releases an energy E comparable 
to that observed in 7-rays, E ~ 10^^ erg, over a time T < lOOs. The large energy density in the 
source results in an optically thick plasma that expands and accelerates to relativistic velocity. 
After an initial acceleration phase, the fireball energy is converted to proton kinetic energy. A 
cold shell of thickness cT is formed and continues to expand with time independent Lorentz 
factor 7 = ~ 300. The GRB is produced once the kinetic energy is dissipated at large radius, 
r > lO^^cm, and radiated as 7-rays through synchrotron and possibly inverse-Compton emission 
of shock accelerated electrons. 



Dissipation may occur due to internal collisions within the shell ( paczyhski Sz Xu 1994 , Rees 
& Meszaros 1994), or due to collision with the inter-stellar medium (ISM) ( [Rees fc Meszaros 
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1992). ISM collision is expected in both cases, since following internal collisions, which convert 
part of the energy to radiation and result from variations in 7 across the expanding shell, the shell 
rapidly cools and continues to expand with approximately uniform Lorentz factor 7 = ??. At first, 
the ISM influence on the shell is negligible. As the cold shell propagates with 7 = t/, it drives a 
shock into the ISM. This shock propagates with 7^ = 2^/^7, and behind it the (rest frame) number 
and energy densities are n' = A'yn and e' = Aj'^nnipC^ respectively, where n is the ISM number 
density ahead of the shock. The width of the shock heated ISM shell is r/47^, where r is the shell 
radius. It is usually assumed that the shell is influenced by the ISM once the heated ISM energy, 
47rr^(r/47^)7^e', is comparable to E. This occurs at a radius ve = {E / A.-Krj^nmpC^Y^^ . However, 
Sari & Piran (1995) pointed out that the ISM may influence the shell at an earlier time, as it 
may drive a relativistic shock wave into the shell at r < te- When such a shock forms the shell 
is no longer "cold". A relativistic backward shock is formed once the energy density e' implies a 
relativistic energy density in the shell, e'/n^rripC^ ~ 1, where n'^ is the cold shell number density. 
This occurs at 

ri? = 2xl0i5(-4^)'^%m, (1) 

where E = lO^^E'sierg, rj = 300t/3oo, n = nicm~^, T = IOTiqs. If > r^;, a relativistic backward 
shock is not formed, and the shell loses its energy to the heated ISM at r^;. If < r^;, most of 
the shell kinetic energy is converted to thermal energy once the backward shock crosses the shell. 
This occurs at (ISari fc Piran 1995| ) 



IQie _±LJ}L cm. (2) 



KAimmpC I \ ni 

At this time the heated shell and ISM move with Lorentz factor 



The ratio of to ve is rn/rE = ^■2hE^[ /T^q %oo- Since -q > 300 is required in order 
for the pair production optical depth to be sufficiently small (e.g. Baring 1993), rji < ve for 
T > Is. In this case, most of the kinetic energy is converted at r = to thermal energy of the 
heated ISM and shell (it can be shown that at this point the shell thickness is similar to that 
of the heated ISM region, rc/47^, and therefore the ISM and heated shell carry similar energy). 
Since the conversion of kinetic energy to thermal energy occurs (as seen by the observer) over a 
time scale of rc/27^ = T, the duration of the GRB is T regardless of whether it is produced by 
internal collisions or by ISM collision. If T < Is and rji > te, only the dissipation at te produces 
the GRB. In this case, the burst duration is VE/ffc. However, it is straight forward to show that 
if we require the duration of the burst to be T^^u, then and the shell Lorentz factor 77 are given 
by the above equations for Vc and 7c with T replaced with T^u- Thus, the conditions at the ISM 
heated shock at the end of the main burst, which determine the following after-glow, are always 
given by (|2|) and (^), with T interpreted as the main GRB duration (the duration over which most 
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of the 7-ray energy is observed). Note that the dependence of rc and 7c on parameters is very 
weak, and that they are mainly determined by E (and the observable T). 

At r > Tc, the shell decelerates and a relativistic shock continues to propagate into the ISM. 



The electrons accelerated during the initial collision at r = rc rapidly lose their energy (Meszaroi 



fc Rees 1997] ). However the shock driven into the ISM continuously produces relativistic electrons 



that may produce the delayed radiation observed on time scales of days to months. As before, the 
shell expanding with Lorentz factor 7 drives a shock with 7^ = 2^/^7 into the ISM. The heated 
ISM energy is 47rr^(r/47^)7'^e' = E/2. This determines the evolution of 7, 

7 = 7c(r/re)-'/'. (4) 

Light emitted from the shell at radius r is observed over duration t = r/2'y'^c = T{r/rc)^. This 
completely determines the evolution of the shell of baryons. In order to calculate the synchrotron 
emission from the heated ISM shell we need to determine the magnetic field and electron energy. 
We assume that the magnetic field (in the shell rest frame) is a fraction of the equipartition 
value, B' = ^^(STre')^/^, and that the electrons carry a fraction of the energy. Since the 
Lorentz factor associated with the thermal motion of protons in the shell rest frame is 7, this 
implies that the Lorentz factor of the random motion of a typical electron in the shell rest frame 
is 7em = '^e7?Tip/me. With these assumptions, and using eqs. (§), (P) and (|^, we find that the 
observed frequency of synchrotron emission from typical electrons, Vm — 77em^-^V'^eC, varies 
with time as 

l^m = IQl'^^^^B^^sf iZaf Hz, (5) 

where t = ltdaydays. This is the frequency where synchrotron emission peaks. Emission at higher 
frequencies is produced by electrons with energy higher than typical. The ratio of synchrotron 
cooling time, tg = Gnrriec/aT^eB''^ , to rest frame expansion time, td = r/'jc, is 

ts/t, = OmC'CB'E^i^'Cyn^'^'. (6) 

From eq. (^, we expect the synchrotron cooling time to be larger than the dynamical 
time for long delay after-glow. As we show in sec. 3, for GRB970228 the parameters are 
such that for most of the after-glow ts > td- In this case, the flux at can be obtained 
in the following way. The number of radiating electrons at a given time is 47rr'^n/3, each 
producing synchrotron photons at a rate (in the shell rest frame) Ps = 'jem'meC^ /tsihi^m/^) (the 
frequency at the rest frame is smaller by a factor 7). Since Vm oc r~^, emission over a range 
Aujn occurs over radius range Ar = Afm,r/6i/m, which corresponds to a time in the rest frame 
At = Ar/'jc. Thus, the energy (in the observer frame) emitted over the frequency range Af^ is 
AE = (27r/9)r^n(7em"ieC^/isc)Afm/i^m- This energy is observed over a time r/2j'^c. Thus, the 
observed intensity at the frequency z/^ is 

^ ^E ^ 2 CTTrrieC ^ jsE 
Av^r/2^^c 27V2^m^/2e^"47rd2' 
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where d is the source distance. If the electron distribution has a high energy tail, dNe/d'je oc 7^ ^' 
for 7e > 7em! then for u > Um (defining a = {p — l)/2): 

F, = F,^[z./z.™(t)]-". (8) 

Note that eq. (^) implies that at early times ts < td (unless the efficiency of converting kinetic 
energy to radiation is very low). We show in sec. 3 that for GRB970228 the parameters are such 
that indeed ts <C at early times. In fact, it has been shown by [Sari, Narayan fc Piran (1996) 
that ts <C td is typically required for the GRB emission. Thus, although p is likely to be similar 
for the GRB and for the after-glow, since in both cases high energy electrons are assumed to be 
produced by shock acceleration, the spectral index a for the GRB is expected to be larger by 1/2. 

Let us now consider inverse-Compton emission from the accelerated electrons. The ratio of 
inverse-Compton to synchrotron emission, Pfc/Ps, is given by Pjc/Ps = Us/Ub, where Us and 
Ub are the synchrotron photons and magnetic field energy density. For t^ < ts, and when emission 
is dominated by the synchrotron process, the total energy in synchrotron photons is a fraction 
td/tg of the electron energy, and the energy density Us is a fraction td/4:ts of the electron energy 
density (the factor of 4 is due to the fact that the rest frame thickness of the shell is r/47, while 
the radiation occupies a shell of thickness r /j) . Using (^) we then obtain 

Pic/Ps = 9ey/'Ei{\^^\ (9) 

The inverse-Compton frequency is a factor 7g larger than the synchrotron frequency, 

= 7 X lO'^ee^Bn^'/'EllS-^/'m. (10) 

Eq. (|9|) implies that > 0.3 is required in order for inverse- Compton emission to be important at 
~ Iday. For ^ = 0.3, however, implies that in order to obtain emission at optical wavelength 
^B < 10~^ is required. It is easy to see from eq. that this would require E to be many orders of 
magnitude larger than the observed 7-ray energy in order for an optical flux to be detectable (see 
also sec. 3 below). Since the observed 7-ray energy is already challenging for most sources, this 
is not likely. Although inverse-Compton emission can not be important when optical emission is 
seen, it will dominate at early times (unless ^ 1), shifting the predicted synchrotron after-glow 
emission to energy well above X-ray/optical. This is discussed in more detail below. 



3. GRB970228 and GRB970402 

Let us first discuss GRB970228. Eq. (^ shows, that the time and frequency dependence 
of is determined by p. For GRB970228, X-ray flux of (2.8 ± 0.4) x lO'^Serg/s cm^ was 
observed lOhr after the burst at the 2-lOkeV band ( [Costa et al. 1997a ), corresponding to 



= {2± 0.3) X lO^^'^erg/s Hz cm^ at = lO^'^Hz, and a flux Fj, = 1.2 x W'^^erg/s Hz cm^ at 
= 3.4 X lO^^Hz {mi = 20.6) was detected after 21hr ( proof et al. 1997[ ). Using these values in 
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we obtain a = 0.63, or p = 2.3. Since the frequency and flux range are large, a is determined 
with good accuracy: The accuracy of the optical fluxes is not quoted; Assuming an accuracy of 
20%, the (la) error in a is ~ 0.05. At 21hr delay, there is a simultaneous detection my = 21.3, 
corresponding to Fi, = 1.1 x 10~^^erg/s Hz cm^ at v = 5.5 x lO^^Hz. This is consistent with 
a = 0.63. However, since the fluxes are probably accurate to only ~ 10%, and since the frequency 
difference is small, this is not a strong test for the model. 

More stringent tests are obtained by observations made at larger delays. With a = 0.63, 
the flux at a given frequency should drop with time as t^^. Following the initial detections, 
observations in both X-ray (2-lOkeV) and optical (V, I) bands were carried out at a time delay 
larger by a factor of 10 compared to the detection delay (87hr for X-ray, 9days for optical). 
Although no flux is quoted for the second X-ray observation, it is mentioned that the flux dropped 
by a factor of nearly 20. Optical observations yielded limits of my > 23.6 and mj > 22.2 at 
9day, implying a drop in flux by a factor larger than 8 an 4 respectively. These observations 
are consistent with the model. HST detected the source in the optical band at 38day delay 
with niv = 26.0 it 0.3 and mj = 24.6 it 0.3 ( pahu et al. 199"7| ). Using the fluxes measured at 
21hr delay, and the scaling, the fluxes predicted by the model at 38day are my = 25.4 and 
mj = 24.7, in remarkable agreement with the observed fluxes (Note that a decline would 

predict my = 27.4, m/ = 26.7). 

The 7-ray fluence of the main GRB is 2 x 10^ erg/cm , emitted over ~ 4s duration with a 
spectral peak at 100 — 150keV (palmer et al. 19"97|) [BeppoSax detected 2 additional sub bursts 
and some lower flux level extending over 80s ( |Costa et al. 1997a| ); however, the peak fluxes of 
subsequent bursts are lower by a factor of 10, and the total emission during the 80s is at most 
twice that of the main burst ( Palmer et al. 199^ )]. This fluence determines the value of f-^E/ATrd? , 



where is the fraction of energy converted to 7-rays {E = 10^^ erg for f^ = l and d = 2Gpc, i.e. 
z ~ 0.5). Using this value in we have 

FJ^ = 1.7 X 10-27^/-ieijerg/s Hz cml (11) 

The value of F^^^ from X-ray and optical detections can be derived from (|8|) if I'mit) is known for 
some t. This can be determined by frequent observations of at a fixed frequency or a by obtaining 
a wide spectrum at fixed time. Such observations are not available for GRB970228. However, we 
can place limits. The optical detections imply that the time tj at which the synchrotron peak is 
at i^rn = 3.4 X lO^'' is t/ < 21hr. Using this in (0) we have 

Fj^ = 1.2 X 10-28 t^i^erg/s Hz cm^. (12) 

where tj = 21t7-^2ihr. This is in remarkable agreement with the 7-ray inferred flux. 
Comparing ( [Tl] ) and (|l^) and using (|5|) we find 
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where d = 2(i2Gpc. Thus, the magnetic field should be close to equipartition, unless the conversion 
of energy to 7-rays is very inefficient. The value of depends mainly on tj. Both and would 
be strongly constrained by observationally determining tj. Using (|l^) in (^) we find that ts > td 
for t > 10^ f^n~^^'^El^(^t J 2i^s. This justifies our assumption that for most of the after-glow we can 
assume tg > t^, and that tg < td at earlier times, with the implications we mentioned for a and p. 

An argument similar to that following ( p!0| ) for the X-ray detection at lOhr, shows that inverse- 
Compton emission can not be important at this time. Using @ this implies < ^.'iE^^^^n-^^ . 
For a given value of the time up to which inverse- Compton dominates the emission is 

t/C = 10niE5i(Ce/0.3)^hr. (14) 

[Eq. dll) is valid for > Q.lhE'^'^^ n~[^'^^ f^'^^ , since we have implicitly assumed that tg > td 
for t > tjc in our calculation, which does not hold for lower values.] For t < tjc, emission 
is dominated by the inverse-Compton process, which produces very high energy photons and 
suppresses the X-ray/optical after-glow. 

For GRB970402 only X-ray after-glow has (so far) been reported ( piro et al. 1997| ). It is 
therefore impossible to constrain the relevant parameters for this burst. However, it is remarkable 
that the detected after-glow is consistent with the model described above. The 7-ray peak flux 
of GRB970402 is 0.08 times that of GRB970228. Assuming both bursts have similar intrinsic 
properties {E, p, and ^b), this implies that the X-ray flux of this burst measured at llhr delay 
(the flux was integrated over 8 to 15hr, and for decline the average flux corresponds to the flux 
at llhr) should be 0.07 times the GRB970228 flux at lOhr, i.e. 2 x 10~^^erg/cm'^s, in agreement 
with the measured flux, 1.5 it 0.5 x lO^^'^erg/cm^s. The non-detection of X-ray emission at 32hr 
delay implies that the flux decreased by more than a factor of 3 compared to the llhr delay, in 
agreement with the model prediction. 



4. Discussion 

We have shown that GRB after-glow observations are naturally explained by the fireball 
model: (i) The functional dependence of the flux on time and frequency is determined only by the 
spectral index p of shock accelerated electrons. GRB970228 after-glow observations are consistent 
with the predicted scaling (^), and require p = 2.3 it 0.1, well within the range expected for 
shock acceleration; (ii) The absolute value of the observe flux (|l^) is consistent with that inferred 
through the model from the 7-ray fluence (0). Apart from determining p, the observed value of 
the flux and the observed normalization of the Um oc relation provide, through eqs. (^) and 
(|7|), two constraints on model parameters {£,e-,^B, E}. These constraints are given for GRB970228 



in eq. (13), where the fireball energy E is replaced with the fraction of energy converted to the 
observed 7-ray fluence, n is the ISM density, and d is the burst distance. The parameter tj that 
appears in (O) is the time at which Um is in the I band, which can only be constrained to tj < 21hr 



for this burst. This parameter may be better constrained by frequent optical observations at 
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shorter delays, or by obtaining a wide spectrum at a single time. This would allow to put strong 
constraints on the least well known model parameters, .^e and (,b- Current observations indicate 
that both can not be much smaller than unity. Note that although the dependence of the .^g 
constraint on d is weak, local GRBs with d ~ lOOkpc instead of d ~ IGpc would imply ^ 1) 
which is of course impossible. Thus, the observed after-glow can not be explained by the model 
considered here for local GRBs. 

Using (^), the derived value of a and 21hr V band observation, the future behavior of 
GRB970228 should follow 

m = mo + 2.51og(t/to), (15) 

where tuq = 25.1, 24.2, 23.7, 22.3 in V, I, J, K bands respectively with Iq = 30day {t measured 
from the GRB). For GRB970402 there is only one X-ray measurement, and it is impossible 
to determine whether the parameters {^e,CB,p} are the same as for GRB970228. If, however, 
these parameters do not vary much between bursts, (|l^) should also describe GRB970402 with 
mo = 24.2, 23.2, 22.8, 21.3 in V, I, J, K bands with to = Iday. 

For the parameters required for GRB970228 the emission is dominated by the inverse-Compton 
process, which produces very high energy photons and suppresses the X-ray/optical after-glow, 
for t < tic = 10(^e/0.3)^hr (for ^ 0.1). The observations therefore imply ^ 0.3. The strong 
dependence of tjc on implies that tic may vary widely from burst to burst. On the other-hand, 
measuring tic by observing the appearance (absence) of X-ray/optical after-glow after (before) 
tic would provide a strong constraint on For ~ 0.2, the inverse-Compton scattered photons 
have a characteristic energy of ~ IGeV during the first ~ 2hr after the burst, therefore providing 
a natural explanation for the delayed GeV emission observed in several strong bursts ([Hurley e\ 



al. 1994)) . X-ray and optical after-glow observations therefore seem to support the suggestion of 



Meszaros &: Rees (1994)| that the delayed GeV emission is due to inverse-Compton scattering of 



ISM electrons heated by the shell collision with the ISM. Note, however, that in Meszaros &: Rees 



(1994) the frequency and duration of inverse-Compton emission is strongly dependent not only on 

but also on the initial fraction of fireball energy carried by rest mass, vic oc rf and tic oc f]'^- 
As we have shown here, the after-glow is independent of rj [cf. eq. (|5|)-(p!o|)]. 

An underlying implicit assumption of the analysis presented in sec. 2, is that ^g, and p are 
constant during the shell deceleration period. Since the ISM shock is relativistic up to a month 
after the GRB event, 7^ = 2{E^i/niY/^{t/mday)-'^/^ [cf. eq. (|), (|), (g)], it is not unreasonable 
that the parameters p, ^g and ^b, which are determined by the physical processes at the shock, 
would remain constant during this time. If these parameters are time dependent, the predicted 
scaling (^ would not hold. The good agreement of GRB970228 after-glow observations with 
indicate that this is not the case. However, it should be noted that at delays much larger than 
months, the shell is no longer relativistic and the predicted scaling laws may no longer hold. 
Another implicit assumption that was made is that the expanding shell is spherical. The results 
are valid, however, also for a jet of opening angle 9, as long as ^ > I/7. At times for which 7 < 1/^ 
the flux would decrease much faster than predicted for a spherical shell. Once again, GRB970228 
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after-glow observations imply that opening angle is not small. It should be emphasized, that since 
dl^ ) imply that both and are not much smaller than 1, a significant decrease with time in 
these parameters, or similarly a decrease in 7 below 1/9, would both result in fluxes much lower 
than observed. 

I thank J. N. Bahcall and P. Kumar for helpful comments. This research was partially 
supported by a W. M. Keck Foundation grant and NSF grant PHY 95-13835. 

REFERENCES 

Baring, M. G. 1993, ApJ, 418, 391 
Costa, E. et al., 1997a, lAU Circular 6572 
Costa, E. et al., 1997b, lAU Circular 6576 
Feroci, M. et al., 1997, lAU Circular 6610 
Piro, L. et al., 1997, lAU Circular 6617 
Groot, P. J. et al., 1997, lAU Circular 6584 
Hurley, K. et al. 1994, Nature, 372, 652 
Meegan, C. A. et al 1992, Nature, 355, 143 

Meszaros, P. 1995, Proc. 17th Texas Conf. Relativistic Astrophysics (NY Acad. Sci., NY) 759, 440 

Meszaros, P., & Rees, M., 1994, MNRAS, 269, 41P 

Meszaros, P., & Rees, M. 1997, ApJ, 476, 232 

Palmer, D., et al., 1997, lAU Circular 6577 

Paczynski, B. 1992, Nature, 355, 521 

Paczyiiski, B., & Rhoads J. 1993, ApJ, 418, L5 

Paczynski, B. k Xu, G. 1994, ApJ, 427, 708 

Piran, T. 1994, in Gamma-ray Bursts, eds. G. Fishman et al. (AIP 307, NY 1994) 
Rees, M., & Meszaros,P. MNRAS, 258, 41P 
Rees, M., & Meszaros, P. 1994, ApJ, 430, L93 
Sahu, K. et al., 1997, lAU Circular 6619 



- 10- 

Sari, R., & Piran, T, 1995, ApJ, 455, L143 

Sari, R., Narayan, R.,&; Piran, T, 1996, ApJ, 473, 204 



This preprint was prepared with the AAS lAT^^X macros v4.0. 



